Stable gapless superconductivity at strong coupling 
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We study cross-flavor Cooper pairing in a relativistic system of two fermion species with mis- 
matched Fermi surfaces. We find that there exist gapless phases which are characterized by either 
one or two gapless nodes in the energy spectra of their quasiparticles. An analysis of the current- 
current correlator reveals that, at strong coupling, both of these gapless phases can be free of 
magnetic instabilities and thus are stable. This is in contrast to the weak-coupling case where there 
are always two gapless nodes and the phase becomes magnetically unstable. 

PACS numbers: 



I. INTRODUCTION 

In recent years, the interest in degenerate fermionic 
systems has considerably increased. In part, this was 
driven by a substantial progress in experimental studies 
of trapped cold gases of fermionic atoms Q. By mak- 
ing use of various techniques, it has become possible to 
prepare atomic systems of different composition, temper- 
ature, density, and coupling strength of the interaction. 
Because of such a flexibility, the basic knowledge gained 
in these studies is likely to be of immense value also out- 
side the realm of atomic physics. For example, the knowl- 
edge of the ground state of an asymmetric mixture of two 
atomic species can shed light on the physical properties 
of strongly interacting dense quark matter that may exist 
in stars. 

It is conjectured that the baryon density in the central 
regions of compact (neutron) stars is sufficiently high for 
crushing nucleons (and strange baryons if there are any) 
into deconfined quark matter. The ground state of such 
matter is expected to be a color superconductor^, Q . 
(For reviews on color superconductivity, see Ref. jj|.) At 
high density the energetically prcferredphase is the so- 
called color-fiavor-locked (CFL) phase j5j]. At densities 
of relevance for compact stars, however, the situation is 
much less clear. 

The difficulties in predicting the ground state of dense 
quark matter are related to the fact that the condi- 
tions of charge neutrality and [3 equilibrium in a macro- 
scopic bulk of matter such as the core of a compact 
star have a disruptive effect on quark Cooper pairing 
0, 0- For example, as emphasized in Ref. en- 
forcing charge neutrality can result in unconventional 
forms of superconductivity, such as the gapless 2-flavor 
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color-superconducting (g2SC) phase. A 3-flavor version, 
the so-called gapless color- flavor-locked phase (gCFL), 
was also proposed 0. While it was argued that both 
types of ga ple ss phases are (chromo-) magnetically unsta- 
ble [nilill[T2| (for the non-relativistic case, see Ref. [f3|). 
there exist convincing arguments that unconventional 
Cooper pairing in one form or another is unavoidable 
(l4| . By taking into account the general observation of 
Refs. [l5L ITfjl \v\, HH . Il9j that the role of gapless phases 
diminishes with increasing coupling strength, one may 
naively conclude that gapless phases do not exist in the 
regime of strong coupling. In this paper, we show that 
this conclusion is premature: not only do such phases 
exist, but they can even be magnetically stable. 

Since the QCD coupling becomes strong at low densi- 
ties, one may conjecture that quark matter undergoes 
Bose-Einstein condensation (BEC) of diquarks rather 
than forming the usual Coo per pairs of the Bardeen- 
Coopcr-Schrieffer (BCS) type jMwl HHH. If this is in- 
deed the case, one may observe well-pronounced diquark- 
pair fluctuations in the vicinity of the critical tempera- 
ture [2 ^. J25I I2f3 | and/or the formation of a pseudogap 
phase |2j. Therefore, it is of interest to study Cooper 
pairing in the strongly coupled regime in more detail. 

In this paper, we address the issue of the existence of 
stable gapless phases in a strongly coupled system of two 
species of massive fermions. (For a related study, based 
on an effective low-energy description, see also Ref. |27j). 
To keep the discussion as simple and as general as possi- 
ble, we consider a model with a local interaction that 
describes two fermion species with equal masses, but 
with non-equal chemical potentials. We find that gapless 
phases with either one or two effective Fermi surfaces can 
exist at strong coupling. These phases are stable in the 
sense that they are free of (chromo-)magnetic instabili- 
ties [iol | . It is expected that many results of our analysis 
should remain qualitatively similar also in the more com- 
plicated case of non-equal masses. 

The paper is organized as follows. In the next section, 
we briefly introduce the model and set up the notation. 
In Sec. IIIII we present the zero-temperature phase dia- 
gram in the plane of the average chemical potential and 
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the coupling constant for several values of the mismatch 
between the fermion chemical potentials. This diagram, 
while obtained in the mean-field (MF) approximation 
and, thus, not completely reliable at strong coupling, re- 
veals a very interesting feature: it has regions of g apless 
phases that are free of the Sarma instability [23 with- 
out imposing the neutrality condition. Note that this is 
drastically different from the situation at weak coupling, 
where the absence of such an instability is mainly due to 
charge neutrality [1, @ , or other types of constraints on 
the system [2||. In Sec. IIVI we calculate the screening 
masses, defined by the long-wavelength limit of the static 
particle-number current-current correlator, and reveal a 
region of parameters (generally at strong coupling) for 
which gapless superconductivity is stable. The discus- 
sion of the results and conclusions are given in Sec. El 

II. MODEL AND QUASIPARTICLE SPECTRUM 

Let us start by introducing the Lagrangian density of 
the model, 

£ = ip(i ? - m + /i7°)V + £i, (1) 

where ip denotes the Dirac field which has two inter- 
nal degrees of freedom, called "flavors" in the following. 
In general, the mass is a diagonal matrix of the form 
m = diag(mi, TO2). For simplicity, however, we restrict 
ourselves to the case of equal fermion masses, i.e., we use 
mi = ni2 = m in this paper. No constraints on the values 
of the chemical potentials of the two flavors of fcrmions 
are imposed, i.e., jl = diag(/ii, ^2) where /xi and \i2 need 
not be equal. 

In order to study superfluidity/superconductivity that 
results from Cooper pairing of different flavors of fcr- 
mions, we introduce the following local interaction term 
to the Lagrangian density, 

d = G(^7 5 ctiC^ t )(^ t z7 5( tiC^), (2) 

where G is the coupling constant, o\ is the (symmet- 
ric) Pauli matrix in flavor space, and C — 17072 is the 
charge conjugation matrix. This term describes a cross- 
flavor attractive interaction between fermions that can 
drive the formation of spin-zero (and, therefore, totally 
antisymmetric) Cooper pairs at weak coupling, or even 
cause the appearance of localized bound states at strong 
coupling. At sufficiently low temperature, these bosonic 
states should form a condensate in the ground state. The 
explicit structure of the condensate is given by the fol- 
lowing expression, 

A = 2G{ip T Ci^aiTp) . (3) 

In the MF approximation, the value of A is determined 
by the minimization of the effective potential 

^ a ) = ^-|S/^MS-™ (4) 



where S 1 (K) is the inverse fermion propagator in 
Nambu-Gor'kov space, 

[S(K)]-i=( * + A7°-™ -ia^A \ 

L v /J V — tai'-fsA _/t — fjrf u —ml K ' 

The corresponding Nambu-Gor'kov spinor is defined by 
\E fT = (ip 7 ' ,4'c) with ipe = Cijr being the charge- 
conjugate spinor. 

The (eight) poles of the determinant det [S(.K')] deter- 
mine the dispersion relations of (eight) quasiparticles. In 
the case of equal fermion masses, these can be given ex- 
plicitly in analytical form, 

fco = ± (e± ± Sy) , (6) 

where all eight sign combinations are possible. In the 
above equation, we used the following notation: 

e± = ^{E k ±pf + A^ (7) 
E k = Vk 2 +m 2 , (8) 



5n = - (10) 

(without loss of generality, we assume that p, > and 
5[i > 0). 

Let us first clarify the structure of the quasiparticlc 
excitation spectrum keeping A/(5/i and fi/m as free pa- 
rameters, i.e., without actually solving a gap equation. 
It is easy to see that, if A/5fi > 1, there are no gapless 
excitations irrespective of the value of jllrn. This is not 
always the case when A/5/i < 1, cf. Ref. In this case, 
gapless modes may exist around effective Fermi surfaces 
at momenta 



k± = \j {p. ± y/sj? - A 2 ) " - m 2 . (II) 

Of course, k± has to be real, otherwise the correspond- 
ing effective Fermi surface does not exist. If fi/m is suf- 
ficiently large, both k + and fc_ are real, and there are 
two effective Fermi surfaces. For A ^ 0, we call the 
corresponding gapless superconducting phase the gSC(2) 
phase. Such gapless phases were discussed in the context 

of dense quark matter @,H[H[ll[ll[ll[lil. A tyP 1 ^ 1 
excitation spectrum is depicted in Fig. ^a). For smaller 
values of p./m, k_ has a non-zero imaginary part while 
fc+ is real, and there is only one effective Fermi surface, 
cf. Fig.^b). Wc call the corresponding gapless supercon- 
ducting phase the gSC(l) phase. Finally, for even smaller 
values of /2/to, both fc_ and k + have non-zero imaginary 
parts, and there are no gapless excitations in the spec- 
trum; we are in a regular, gapped superconducting (SC) 
phase, cf. Fig. ^c), just as in the case A/8fi > 1, cf. 

Fig.md). 

In Fig. [3 wc show where the three different types of 
superconducting phases occur in the plane of A/S/j, and 
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FIG. 1: Typical quasiparticle excitation spectra in gapless 
and regular superconducting phases. The values of \xjm and 
A/Sfj, are specified inside the panels. In all four cases S/j,/m = 
0.5. Excitation spectra of antiparticles are not shown. 



p/m, for a fixed value of dfi/m = 0.5. The boundary of 
the region gSC(l) can be derived from the requirement 
that k- = 0, i.e., the region exists for values of A/S/j, and 
p/m satisfying the condition 



/! — m 
8fi 



< 1. 



(12) 



The boundary between the gSC(2) and SC regions is 
given by A/Sfi = 1 for p/m > 1. The three regions 
SC, gSC(l) and gSC(2) meet at the "splitting" point S, 
cf. Ref. Wi- 



lli. PHASE DIAGRAM 

In this section, we explore the phase diagram of the 
model at hand in the plane of the average chemical po- 
tential p and the coupling constant G. We shall restrict 
our study to the zero-temperature case when the prob- 
lems with the stability of gapless phases are expected to 
be most prominent. Then, the effective potential reads 



A 2 

- 4G 



k 2 dk 
2tt 2 



6n\ 



■5fi\ 



(13) 



where A is a momentum cut-off. 

In Fig. [3J we show the phase diagram for three differ- 
ent values of the fermion masses, m/A = 0.1, 0.2 and 0.3, 



to 




u/m 



FIG. 2: The phase diagram in the plane of A/Sfi and p/m 
plotted for a fixed value of Sfi/m = 0.5. The region SC repre- 
sents a superconducting phase without gapless modes, while 
gSC(l) and gSC(2) represent phases with one and two effec- 
tive Fermi surfaces, respectively. The three regions merge at 
the splitting point "S", cf. Ref. l27ll . For a given m/A, the 
region to the right of the dashed line (from top to bottom, 
m/A = 0.3, 0.2, and 0.1) is not accessible in the MF analysis 
of Fig. |3| see text for a detailed explanation. For m/A = 0.2, 
the "screening mass" uim is imaginary in the shaded area. 



and a fixed value of Sfi/m = 0.5. The coupling constant 
is normalized by Go = 47r/A 2 . Note that, for m = 
and for G > Go, the vacuum is unstable with respect to 
BEC of diquarks. The bold and thin solid lines repre- 
sent first- and second-order transitions, respectively For 
each choice of mass, the part of the diagram above the 
transition lines corresponds to the fully gapped super- 
conducting phase. The small regions bounded by the 
solid and dotted lines in the middle of the phase diagram 
correspond to gapless phases. Below the solid lines, the 
system is in the normal phase. The normal phase to the 
left of the thin dashed vertical lines at p — m — Sfi cor- 
responds to the vacuum. 

As one can see from Fig. [21 at weak coupling the 
phase transition between the normal and superconduc- 
ting phase is of first order. This is quite natural when 
there is a fixed mismatch between the chemical poten- 
tials of pairing fcrmions. With increasing the coupling 
constant G, the critical chemical potential p becomes 
smaller and the lines of first-order phase transitions ter- 
minate at endpoints. It is worth mentioning that the 
endpoints lie completely inside the region of the super- 
conducting phase when m/A = 0.1 and 0.2. In the case 
m/A = 0.1, one can see this more clearly from the in- 
sertion in the lower left corner of Fig. [31 In the case 
m/A = 0.3, however, the corresponding endpoint lies on 
the phase boundary between normal and superconduc- 
ting matter, or so close to it that our numerical resolution 
is not sufficient to make a distinction. 

Across the boundary of a second-order phase transi- 
tion, the gap A must be continuous. This means that 
the gap assumes arbitrarily small values just above the 
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FIG. 3: The phase diagram in the plane of G/Go and fi/A 
for three different values of the fermion masses, m/A = 0.1, 
0.2 and 0.3, and a fixed value of Sfi/m — 0.5. The bold 
and thin lines represent first- and second-order transitions, 
respectively. The insertion is an enlarged view of the endpoint 
of the first-order transition line in the case of m/A = 0.1. 



thin solid lines in Fig. [3] When the mismatch between the 
chemical potentials Sfi is non-zero, there inevitably exists 
a region adjacent to the transition line where A < Sfi. 
From Fig. [3 we see that this is a necessary condition to 
have a gapless phase, and this is precisely what we ob- 
serve in Fig. |21 It is, however, not a sufficient condition: 
if fi/A is sufficiently small, we could also have a gapped 
phase, as seen to the left of the gSC(l) region in Fig. [3 
The gapless phases in Fig. correspond to global min- 
ima of the effective potential and, thus, they are free of 
the Sarma instability . This might be surprising since 
we do not impose additional constraints such as neutral- 
ity. An apparent discrepancy between this finding and 
that of Rcf. Q is resolved by noting that the stable gap- 
less phases in Fig. [21 always occur at strong coupling. 

We now want to clarify which regions in Fig. [2 are ac- 
cessible by the MF calculation that gives rise to the phase 
diagram in Fig. |21 We first note that the region below 
the solid lines in Fig. [21 represents the vacuum or the 
normal-conducting phase. Since there A = 0, moving 
along its upper boundary corresponds to moving along 
the horizontal axis in Fig. [3 Across a first-order phase 
transition, the gap is discontinuous. Therefore, certain 
values of A/S/j, are excluded. These can be found by 
computing the values of the gap along both sides of the 
bold solid lines in Fig. [21 For fixed m/A, such a path 
starts at the merging point of the second- and first-order 
transition lines in Fig. [21 continues along the lower side 
of the bold line, and runs around the endpoint, before 
continuing along the upper side of the first-order transi- 
tion line. In Fig. this path is shown as a dashed line. 
The region to the right of this line is not accessible in the 
MF analysis. As a consequence, the only gapless phases 



appearing in the phase diagram in Fig. [21 are those of 
type gSC(l). This excludes, therefore, possible ground 
states that correspond to the gSC(2) region as well as 
the splitting point S. In fact, the statement regarding 
the splitting point can be made rigorous by noting that 
the second derivative of the effective potential d 2 V/dA 2 
is negative at S, meaning that this point cannot be a 
minimum of V(A). Of course, this conclusion may easily 
change if an additional constraint (e.g., such as neutral- 
ity) is imposed on the system. 

Now, it is natural to ask whether the gSC(l) and 
gSC(2) types of gapless phases are subject to the chro- 
momagnetic instability [lfj. 
the next section. 



This is studied in detail in 



IV. STABILITY ANALYSIS 

In this section, we discuss the stability of the gapless 
phases that were introduced in Scc.[n] see Fig. [3 To this 
end, we calculate the fermion-number current-current 
correlator and study when such a correlator points to- 
ward an instability. Note that the use of the fermion- 
number current is not accidental here. When the vacuum 
expectation value in Eq. J2J) is non-zero, the fermion- 
number symmetry is spontaneously broken. 

By definition, the current-current correlator is given 

by 



n""(p) = 



T 



d 3 k 



t^S(K)t v S(K - P) 



2 ^ J (2tt) 3 

(14) 

where T M = diag(i7 M , —17^) is the vertex in Nambu- 
Gor'kov space. Following the approach of Ref. [lj}, we 
consider IP" only in the static (po = 0) and long wave- 
length limit (p — > 0). Then, it is convenient to introduce 
the "screening masses" which arc defined by 



i 2 D = - lim n ou (w = 0,p), 



(15) 



«4r = -o lim >y +PiPj)n y (o; = 0,p). (16) 
z p— >o 

If the fermion-number symmetry is promoted to a gauge 
symmetry, the quantities and tum would describe 
the electric and magnetic screening properties of a su- 
perconductor. 

Our calculations show that m D is positive definite in 
the whole of plane of A/5fi and p/m, see Fig. [2 How- 
ever, the magnetic screening mass mjj could be imagi- 
nary (i.e., m 2 M < 0) in some cases when is nonzero. In 
general, an imaginary result for tum indicates an insta- 
bility with res pect to the formation of inhomogencities 
in the system j3fl l3lj . In the weak-coupling limit, the 
instabilities develop when A < 8/j, [ljj- In this paper, 
we study whether a similar instability also develops at 
strong coupling. 

When A is nonzero, the Nambu-Gor'kov propagator of 
fermions has both diagonal and off-diagonal components, 
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see Eq. |J5J. Their contributions to the current-current 
correlator (|14|> can be considered separately. Then, the 
expression for m 2 M can be given in the following form: 

m 2 M = (m 2 M ) diag + (m 2 M ) oS . (17) 

Each of the two types of contributions in this equa- 
tion can be further subdivided into particle-particle (pp), 



antiparticle-antiparticlc (aa), and particle- antiparticle 
(pa) parts, i.e., 

K/W = (ti&)£2 + (m 2 M ) { Zl + (m^Sfi, (18) 

and a similar representation holds for (m 2 M ) R. The ex- 
plicit expressions for all contributions are given by 



( m 2 N(pp),(aa) 

( m 2 N(pp),(aa) 
\ m M)oB 



(m 2 1 (pa) 

l TO A/Jdiag 



/ 2 \(pa) 
\ m M)oS 



d 3 k k 2 

d 3 k k 2 
d 3 k 

d 3 k 
(2tt) 3 



-_[l-0(-e T + <^)]- 



A 2 

7r[i-0(-*F 

fc 2 



3 ^ 

3-* 



1 



A 2 



A 2 

^)] - — 6 ( e T - 5 t- 1 ) 



e 2 - E+E- 



[1 - 6(-e + + 8n)] - (e+ 
+ Sfi)] - (e+ -> e_) 



(19) 
(20) 



"i^' (21) 



(22) 



where S± = ± <5/i. The upper and lower signs in 
Eqs. (fH3|) and (f^TTf) denote the (pp) and (aa) parts, re- 
spectively. It should also be noted that the vacuum con- 
tribution 1/Ek was subtracted in Eq. (|21|) . 

Several remarks are in order regarding the expressions 
in Eqs. (|T§jl through Firstly, we note that the con- 
tributions due to the terms with the 9- and <5-functions 
in the integrands vanish in the gapped superconducting 
phase (SC). These are nontrivial, however, in the gap- 
less phases when there exists at least one well-defined 
effective Fermi surface, see Eq. (|ll|l . Secondly, the terms 
with the 5-function give negative contributions that di- 
verge as — [(^m) 2 — A 2 ] when A — > 5fi from below. 
This is a reflection of the divergent density of states at 
k = k± (see, e.g., the second paper in Ref. @). Finally, 
we point out a qualitative difference between the inte- 
grands in the (pp) and (aa) expressions and those in the 
(pa) ones when k — > 0. Because of the overall factor 
k 2 /E 2 in the former and the factor (3 — k 2 / E 2 ) in the 
latter, the low-momentum contributions of the (pp) and 
(aa) parts are suppressed, while those of the (pa) parts 
are not. This has important consequences at strong cou- 
pling. 

Our numerical results for function of A/SfJ, for 

three different values of the ratio JX/m (i.e., 0.9, 1.3 and 
2) and a fixed value of the mass, m/A = 0.2, are shown in 
Fig. 0] (Note that the qualitative features of the numer- 
ical results are robust when the parameter m/A changes 
in a relatively wide range.) The solid lines represent the 
complete results for the screening masses squared, m 2 M , 
while the other four lines give the separate (pp) and (pa) 
contributions, defined in Eqs. (|19fl through 1221) . We do 
not show the (aa) contributions in Fig. 0] because they 
are always small numerically. 



In each panel of Fig. 2] the value of m 2 M vanishes in the 
limit A/Sfi — > which corresponds to the normal phase. 
As one can infer from the figure, this is due to the ex- 
act cancellation of the negative (paramagnetic) (wtf)diag 
and positive (diamagnetic) ("i 2 /)^ag contributions. 

As one can see from the top panel of Fig.0J the result 
for m 2 M is positive definite in the case p,/m = 0.9. This is 
not so, however, when the value of the ratio Jxjm is larger 
than 1. In particular, the quantity m 2 M develops a nega- 
tive divergence just below A/Sfi — 1 and stays negative 
for a range of values of A/Sfi, see the middle and bottom 
panels in Fig. 0] It is easy to figure out that the diver- 
gence is caused by the singular behavior of the density of 
states around the effective Fermi surfaces in the gapless 
phases. This singularity affects only the particle-particle 
contributions, i.e., (to 2 /)^^ and (m?MhM ■ When the 
value of p,/m becomes larger than 1 and increases further, 
negative values of m 2 M first appear only near A/6fi = 1 
(see, e.g., the middle panel in Fig. When the ratio 
jS/m gets larger than about 1.8, however, the quantity 
m 2 M becomes negative in the whole range A/Sfj, < 1 (see, 
e.g., the bottom panel in Fig.0J). 

The numerical results for the screening mass can be 
conveniently summarized in Fig. [3 by identifying the re- 
gion in which tum is imaginary (i.e., m 2 M < 0). For the 
given set of parameters, Sfi/m = 0.5 and m/A = 0.2, this 
is marked by the shaded area there. For p,/m > 1.8, the 
gSC(2) type gapless phase is unstable in the whole range 
A/S/j, < 1 where it is defined. This is in ag reement, of 
course, with the results at weak coupling 10] . It is most 
interesting, however, that for smaller values of the ratio 
fi/m, the quantity m 2 M could be positive even in the gap- 
less phases. In particular, m 2 M > in the whole gSC(l) 
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FIG. 4: The dependence of m\j (solid line) as well as the 
(pp) and (pa) contributions, denned in Eqs. 119H through 1221 k 
versus A/S/j, for three different values of Ji/m and a fixed value 
of m/A = 0.2. 



region as well as in a part of the gSC(2) region. Note 
that the gapless phases found in the diagram in Fig. 
always correspond to the stable region. 

To complete the analysis of the stability of gapless 
phases at strong coupling, we note that a general require- 
ment is that the eigenvalues of the susceptibility matrix, 
—d 2 V/d[iidiij, are non-negative j^. As a criterion for 
stability, the susceptibility matrix is meaningful only in 
the ground state defined by that A(jS, 6/j.) which solves 
the gap equation. In a sense, therefore, the susceptibil- 
ity criterion has a less general status than the one of the 
(chromo-) magnetic stability. The latter does not refer to 
a specific form of the gap equation. This difference might 
be very important when additional constraints (e.g., neu- 
trality) are enforced on the system. 

In the MF approximation used in this study, the eigen- 
values of the susceptibility matrix are non-negative in all 
phases of Fig. |3 as well as in the regions to the left of 
the dashed lines in Fig. [21 We argue that this is related 
directly to the fact that the ground state is defined as the 
global minimum of the effective potential. Indeed, one of 
the important necessary conditions for the susceptibility 
criterion to be satisfied reads [331 



2 V 



> 0. 



(23) 



d 2 V 
dA 2 



>0, 



(24) 



d{5^) 2 

In the model at hand, as one can easily check, this is 



which is always satisfied at the global minimum of the 
effective potential. 

Before concluding this section, we would like to briefly 
discuss the properties at the splitting point S in Fig. 
This point describes the situation when the two effec- 
tive Fermi surfaces merge exactly at zero momentum. 
This coincides with the definition in Ref. [27j where the 
properties near S are studied using an effective theory 
approach. From Fig. [21 we see that four qualitatively 
different regions merge at point S, i.e., (i) gSC(l), (ii) 
SC, (iii) gSC(2) with a positive value of m 2 M: and (iv) 
gSC(2) with a negative value of rn 2 M . This topology is 
in a gree ment with the qualitative picture presented in 
Ref. although here we use a different approach. It is 
interesting, though, that the ground state defined by the 
splitting point is never realized in the MF approximation 
in our model. This is because the second derivative of 
the thermodynamic potential is negative definite at the 
point S, which excludes the possibility of having a mini- 
mum there. 



V. DISCUSSION AND SUMMARY 

In this paper, we explored the superconducting (su- 
pcrfluid) phases in a relativistic model with two fermion 
species having mismatched Fermi surfaces. We found 
that, in general, there could exist gapless phases with 
either one or two effective Fermi surfaces. In the MF ap- 
proximation, however, only the gapless phase gSC(l) is 
realized as a ground state, and only at strong coupling. 
We also calculated the fermion-number current-current 
correlator in the static and long wave-length limit. The 
results show that the gapless phases at strong coupling 
are free of the (chromo-)magnetic instability. This is in 
contrast to the situation in neutral quark matter at weak 
coupling [Tcj . 

It is interesting to note that there could exist different 
regimes of BEC in the model at hand. In general, bound 
bosonic states are formed when Ji/m < 1 [32^. As seen 
from Fig. [3 this is satisfied in parts of the SC and gSC(l) 
regions, giving rise to gapped and gapless phases, respec- 
tively. While the number densities of the two species of 
fermions are equal in the former, there is an excess of 
one of the species in the latter. In fact, this excess is an 
order parameter that defines the gapless phases at zero 
temperature 0. The gapless BEC phase is a mixture 
of tightly-bound bosons in the form of a condensate and 
additional unpaired fermions |27l |33( . 

It is worth mentioning that the gapless phase of the 
gSC(l) type could also exist for a range of parameters 
when Ji/m > 1. In this case, stable bosons do not exist 
and, thus, the BEC regime cannot be realized. If the di- 
quark coupling in quark matter is sufficiently strong, this 
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gaplcss phase can potentially be realized as the ground 
state of baryon matter. Our analysis of the one-loop 
fermion contribution to the current-current correlator 
suggests that the gSC(l) phase is stable. In the case of 
quark matter, however, it might be important to check 
if the inclusion of gluon and ghost contributions does 
not change the conclusion. While such contributions are 
negligible in the weak-coupling limit, this may not be the 
case at strong coupling. 

A few words are in order regarding non-relativistic 
models [53, H3, 03 HH UK- I n our analysis, we saw 
that purely relativistic effects due to the antiparticle- 
antiparticle loops were always small near the splitting 
point. This suggests strongly that our results should re- 
main qualitatively the same also in the non-relativistic 
limit that is defined by |/2 — m\/m <C 1, Sfi/m. <C 1 and 
A/m <C 1. Indeed, these conditions define a narrow area 
around the vertical line p./m = 1 in Fig. [21 which includes 
the point S. In fact, this also suggests that the topology 



around the splitting point is the same in both relativistic 
and non-relativistic models. 

Note added. While finishing this paper, we learned 
that a similar study in a non-relativistic model is done 
by E. Gubankova, A. Schmitt, and F. Wilczek |37j . 
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